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^ ■ Abstract 
^ ■ 

, The braneworld modeL in which our Universe is a three-brane embedded in a five-dimensional 

(N . 

bulk, allows the existence of wormholes, without any violation of the energy conditions. A fun- 

! 

^ ■ damental ingredient of traversable wormholes is the violation of the null energy condition (NEC). 

CO . However, in the brane world models, the stress energy tensor confined on the brane, threading 

' the wormhole, satisfies the NEC. In conventional general relativity, wormholes existing before in- 
U ■ 

Q^' flation can be significantly enlarged by the expanding spacetime. We investigate the evolution of 

an infiating wormhole in the brane world scenario, in which the wormhole is supported by the 

i nonlocal brane world effects. As a first step in our study we consider the possibility of embedding 

I a four-dimensional brane world wormhole into a five dimensional bulk. The conditions for the 

• embedding are obtained by studying the junction conditions for the wormhole geometry, as well 

, as the full set of the five dimensional bulk field equations. For the description of the infiation we 

I adopt the chaotic inflation model. We study the dynamics of the brane world wormholes during 

■ the exponential infiation stage, and in the stage of the oscillating scalar field. A particular exact 

> '. 

solution corresponding to a zero redshift wormhole is also obtained. The resulting evolution shows 



, that while the physical and geometrical parameters of a zero redshift wormhole decay naturally, 

a wormhole satisfying some very general initial conditions could turn into a black hole, and exist 
forever. 
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I. INTRODUCTION 

Wormholes are hypothetical connected spacetime , and are primarily useful as 

"gedanken-experiments" and as a theoretician's probe of the foundations of general rela- 
tivity. The wormhole has no event horizon, which allows a two-way passage of energy or 
light signals. It is possible to transform wormholes into time machines for backward time 
travel 3|, 1^ . Wormhole solutions are a specific example in solving the Einstein field equation 
in the reverse direction, namely, one first considers an interesting space-time metric, then 
one finds the matter source, responsible for the respective geometry. In this manner, it 
was found that some of these solutions possess a peculiar property, namely "exotic matter" , 
involving a =tre=s-e„ergy tensor that violates the null energy condition. A number of specific 
solutions with this property have been found ([5] and references therein). These geometries 
also allow closed timelike curves, with the respective causality violations [6|]. However, the 
existence of wormholes without horizon would lead to the violation of the energy conditions 
at the throat of the wormholes [3, S]- This implies that for some observers wormholes would 
require negative energy. 

The conventional manner of finding wormhole solutions is essentially to consider an in- 
teresting space-time metric, and then to derive the stress-energy tensor components. A 
more systematic approach in searching for exact solutions, namely, by assuming spherical 
symmetry, and the existence of a non-static conformal symmetry, was considered in [9]. 

The requirement of negative energy makes the existence of a considerable number of 
wormholes nowadays very unlikely. However, the existence of wormholes in the early Uni- 
verse seems to be more natural, as the quantum mechanical effects would be much stronger. 
Quantum mechanical effects may provide the negative energy, like, for example, via the 
Casimir effect. An inflating wormhole in conventional general relativity was already pro- 
posed and studied by Roman [10]. However, the systematic analysis of the cosmological 
implications of the existence of the wormholes in the early Universe still remains an open 
problem. 

While the negative energy required by a wormhole to exist is the main problem of the 
theory in standard general relativity, the braneworld models of the Universe provide a natu- 
ral way for the existence of the wormholes. In this work, we consider the Randall-Sandrum 
brane world model, which were originally introduced to give an alternative to the compacti- 
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fication of the extra dimensions 
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12| , and to explain the hierarchy problem of the particle 



physics. Standard 4D gravity can be recovered in the low-energy limit of the model, by 
assuming that a 3-brane of positive tension is embedded in a 5D anti-de Sitter bulk. The 
covariant formulation of the brane world models has been developed in [isj], leading to the 
modification of the standard Friedmann equations on the brane. It turns out that the dy- 
namics of the early universe is altered by the quadratic terms in the energy density and 
by the components of the bulk Weyl tensor, which both give a contribution in the energy 
momentum tensor. This implies a modification of the basic equations describing the cos- 

n 

mological and astrophysical dynamics, which has been extensively considered recently |14j |. 
There is a number of works that use the brane world model to solve various cosmological 



problems 
ered in 
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or instance, matter exchange between the brane and the bulk has been consid- 



16| . and this gives a possible scenario of reheating [l7||. Scalar field inflation 



181], or non-scalar field inflation 19|, have also been considered in the brane world models. 

In a brane world Universe wormholes could exist without negative energy. This was 
pointed out by Lobo [20|], who showed that due to the braneworld corrections to the energy- 
momentum tensor a static wormhole in a RSII brane could obey the energy condition ev- 
erywhere. The terms that allows the energy condition of the wormhole to be satisfied in 
the brane world model come from the nonlocal projected Weyl tensor, and they appear as 
a correction terms in the AD Einstein equation. The braneworld correction in energy is 
equivalent to ordinary matter in convention al g eneral relativity 21 1. 

An interesting question, raised first in 22|, is the possibility of wormhole-black hole 
transition. It is quite likely, in view of some recent studies, that a time-dependent equation of 
state had caused the Universe to evolve from an earlier phantom-energy model. In that case 
traversable wormholes could have formed spontaneously. Such wormholes would eventually 
transform into black holes. This would provide a possible explanation for the existence of a 
huge number of black hole candidates, while any evidence for the existence of wormholes is 
entirely lacking. 

It is the purpose of the present paper to consider the general properties and the dynamics 
of an inflating wormhole in the brane world scenario. As a first step in our study we consider 
the problem of the embedding of a braneworld wormhole into a five dimensional bulk. By 
considering a general five dimensional metric, with all metric tensor components function of 
time, radial coordinate and the extra-dimension, we formulate first the junctions conditions 
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that must be satisfied by the wormhole geometry. Then, by considering the full system of 
five-dimensional Einstein field equations we obtain the five-dimensional equations satisfied 
by the wormhole metric tensor components. Particular cases of embedding (static geometry, 
homogeneous energy-momentum tensor and scalar field dominated braneworld models) are 
considered in detail, and it is shown that for each of these cases wormhole models can be 
obtained, at least in principle, in a full five- dimensional setup. As a next step in our study we 
consider the evolution of an inflating (four-dimensional) braneworld wormhole in the early 
Universe. The inflation model we choose is the chaotic inflation model in the RSII braneworld 
model, considered in 23|. Current observations still permit the quadratic potential chaotic 



inflation model in the braneworld. On the contrary, the quadric and higher order exponents 
are excluded in the high energy regime. We study the evolution of a braneworld inflating 
wormhole in a chaotic inflation model, with quadratic potential, i.e. V{(f)) = m^(/)^/2, where 
m is the mass of the inflaton field, and is the scalar field. We assume that during inflation 
a braneworld wormhole that does not violate the energy conditions is created by the scalar 
field, or by the local perturbations of the scalar field. We consider how such a wormhole 
would evolve as the Universe inflates, reheats, and becomes radiation dominant. Finally, we 
consider the conditions under which the transition of the wormhole to a black hole could 
have taken place in the early Universe. 

The present paper is organized as follows. In Section [TTl we consider the general prob- 
lem of the embedding of a braneworld wormhole into the bulk geometry. Particular cases 
of embedding are considered in Section UTTl In Section HV] we review, following {20 1, the 
properties of the braneworld wormholes, and show that they can exist without the viola- 
tion of energy condition. In Section |Vl we introduce the inflating braneworld wormhole 
model. The basic equations describing the braneworld evolution of the inflating wormhole 
are presented in Section IVll A simple solution of the field equations corresponding to a zero 
redshift wormhole, is obtained in Section IVIII The chaotic inflation on the brane is reviewed 
in Section IVIIII The evolution of the inflating wormhole during the exponential inflation 
and oscillating phases is analyzed in Section [IXl We discuss and conclude our results in 
Section |Xl In the present paper we use the natural system of units with G = c = 1 
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II. THE BULK CONFIGURATION OF THE BRANEWORLD WORMHOLES 



The four-dimensional wormhole models we are going to consider rely only on the non- 
local projected Weyl tensor of the bulk. Hence it would be important to investigate what 
kind of bulk configuration can lead to a Weyl-tensor supported four- dimensional wormhole. 
The wormhole on the brane could possibly arise due to a specific form of the Weyl tensor 
of the bulk, like, for example, the one corresponding to a gravitational wave, or it may 
arise because our brane is located in such a way that the homogeneity is broken, so that 
the bulk metric is "abnormally" projected. Of course the specific form of the Weyl tensor 
supporting a four-dimensional wormhole can also be the combination of both previously 
mentioned effects. Therefore we can consider the full gauge freedom of the metric of the 
bulk, such that it projects as a wormhole on the brane, i.e., the metric coefficients qab on 
the bulk could all be non-zero. In this way we may think on the wormhole as an effect of a 
5D gravitational wave. However, as motivated by the fact that an expanding Universe can 
be modeled as a motion of the brane in the bulk, we would like to model the brane wormhole 
by allowing the brane to offset its position. 

In the following we denote by capital letters the coordinates in the bulk, with values 
running from to 4. Greek letters mean coordinates on the brane, while lower case letters 
mean spatial coordinate. We denote by the normal vector to the brane, and by y the 
fifth coordinate. The derivatives with respect to the coordinates are denoted by a comma. 



A. The induced metric and the junction conditions for braneworld wormholes 

In a suitable system of coordinates, the expansion of the Universe can be viewed as the 
motion of the brane through a static bulk 2J, |25|. Suppose a region of the brane slowly 



expands in a non- homogeneous manner. This creates a local inhomogeneity. However, we 
assume that asymptotically the expanding Universe is homogeneous. Due to its motion, the 
brane will now be located into a new position 

y = Yir), (1) 

where we assume that Y has no dependency on t. We assume that the 3- homogeneity of the 
bulk metric is broken locally. Hence the Gaussian coordinate system may not be the normal 
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coordinate system locally. Therefore for the bulk metric we assume the general form 

ds"^ = -M{t, r, yfdt"^ + N{t, r, y^dr^ + P{t, r, y)^ [dfi + sin^ Odcj)'^) + Q(t, r, (2) 

with all the metric coefficients functions of the time and radial distance coordinate, as well as 
of the five-dimensional coordinate y. The metric tensor coefficients satisfy the bulk Einstein 
equation 

(^)G^B = -(^)A(^)^?AB + kl [(5)T^B + 5{y - F(r))r^r1 , (3) 

where *^^^A is the five- dimensional cosmological constant, ^^^Tab is the bulk matter energy- 
momentum tensor, and T^^"'' is the energy-momentum tensor of the matter on the brane, 
respectively. The five- dimensional gravitational coupling constant is denoted by k"^. 

If we use t, r, 6, (p as brane coordinates, the normal vector n can no longer be a constant 
throughout the brane, since the brane is offset. The induced metric on the brane is obtained 
by substituting y = Y(r) in the metric coefficients, and by taking into account the presence 



of Qdy in the bulk, as 



26| 



= -M(t,r,F(r))2dt2+ [iv(t,r,r(r))2 + g(t,r,r(r))%(r)2]rfr2 + 

P(t, r, F(r))2 {de'^ + sin^ ed(f)'^) . (4) 

In particular, this metric will induce an infiating wormhole if the metric coefficients satisfy 
the following conditions 

M{t,r,Y{r)y = e^*^^'*), (5) 

N{t, r, y (r))2 + Q(t, r, F (r))%(r)2 = [ , (6) 

1 — b{r)/r 

P{t,r,Y{r)f = a{tfr\ (7) 

g24>(r,t) -g Qg^iig(j |;]2g redshift function, and h{r) is called the throat function of the worm- 
hole. The throat of a wormhole is defined by 6(ro) = tq. A time dependent throat function 
would lead to divergence of energy density at the throat and therefore we do not consider 
it. The function a{t) is the scale factor describing the expansion of the wormhole. 

By integrating the bulk Einstein equations along the extra dimension it follows that the 
matter content in the brane affects the bulk metric through the iump of the brane extrinsic 

n 

curvature K^y [27[. This is the junction condition of the brane. The extrinsic curvature can 
be calculated according to 

Kab = h^VcUB, (8) 
6 



where 



hAB = Qab - rtAUB, 



(9) 



is the projection tensor which is a tensor in the bulk that acts as metric on the brane. The 
brane is assumed to have the same properties on both sides. Mathematically this represents 
the Z2 symmetry. By imposing the Z2 symmetry we can relate the extrinsic curvature to 
the energy-momentum tensor of the matter on the brane, 



K,,{r,t,R{y+)) - K^,{r,t,R{y-)) = -2K,,{r,t, R{y)) 



T, 



fii/ 



iTh 

3 



(10) 



where the second equality follow from the Z2 symmetry of the brane, and the metric and the 
energy-momentum tensors of the brane can be obtained from the corresponding tensors of 
the bulk by using the coordinate transformation defined by Eq. ([1]) . The extrinsic curvature 
specifies the derivative of the metric along the flow defined by the normal vector. It tells us 
how the metric evolves off the brane in order to match the external space. It also sets the 
boundary conditions for off brane evolution. The components of the normal vector can be 
calculated according to 

ua 



where 



which give 



U = ^A {y-Y{r)) 



Ua 



0, 



.,0,0, 



1 



(11) 



(12) 



(13) 



'1/Q2 + Y^JN^ ^l/g2 + y2/jv2 

The energy - momentum tensor of the matter on the brane is given by a delta function 
5{y — Y{r))T^Q^'^. For the brane energy - momentum tensor T^^^'^ we assume the perfect 
fluid form, 



T 



brane 
AS 



[p(t, r) + p(t, r)] UaUb + p(t, r)hAB, 



(14) 



with p(t, r) and p(t, r) are the energy density and the pressure of the matter on the brane, 
respectively. Ua is the five - dimensional velocity field of the matter fluid on the brane. 
T^B^"^ satisfies the condition T^^^'^n'^ = 0. The four- dimensional energy- momentum tensor 
Tf^u can be obtained from T^^"'' by the relation T,,,u = e^e^TAs, where are the tetrad 
vectors that form an orthogonal basis on the brane 26| . In the following we assume that the 
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fluid is at rest with respect to the comoving observers, so that Ua = [e*, 0, 0, 0, 0]. Therefore 
the junction conditions can be written as 

" NQ^/WTYJQ^ 6 ' ^ ^ 

-Y, {Q,N - N,Q) 
Ktr = ^2 + y2g2 = 0' (16) 

y..Q^iV^ - Y^^QQ.N^] = -f '^,^y,f' , (17) 



^'^e - ^r^/Ar9 , ^^9^9M/9 - IT^ P' (l^J 



Arg(Ar2 + y2g2)i/2 " g 

K^^ = sin^ ^i^e,. (19) 



In order for a brane- world wormhole be embedded into the bulk these junction conditions 
must be satisfied. The first condition for the embedding is obtained from Eq. f lTC]) as 

NAt.r.Yjr)) ^ g,(t,r,y(r)) 

N{r,t,Y{r) Q(r,t,F(r))' ^ ^ 

since we have assumed that Yr 7^ 0. In a brane without radial energy flow, a co-expanding 
wormhole that travels in the y direction imposes the condition that the metric coefficients 
and Q have a similar time behavior on the brane. By substituting Eq. (!20|) into the time 
derivative of Eq. (Q we obtain a separable equation with solutions of the form 

N{t,r,Y{r)) = a{t)Af{r),Q{t,r,Y{r)) = a{t)Q{r), (21) 

where A/'(r) and Q{r) are some arbitrary functions of the radial coordinate only. The position 
of the brane Y{r) on which an inflating wormhole can exist is given by the equation 



YHr) = — 



(22) 



1 — b{r)/r 

The throat of the wormhole can be pictured as a travel upward relative to the brane, i.e., 
y(ro) —7- 00. Due to the Codazzi equation for the extrinsic curvature VbK^ — VaK = 



8 



RbcSa^'^^ junction condition also implies the local conservation of the energy- 
momentum tensor. On a source free bulk we obtain the standard conservation equation 
of the energy-momentum tensor on the brane, 

= 0, (23) 

which gives the following constraints on a and $ with respect to the evolution of the matter 
on the brane, 

p^t{t, r) + 3- [p{t, r) + pit, r)] = 0, (24) 
a 

9,$ [p(t, r) + pit, r)] + drpit, r) = 0. (25) 
In the above equations we have used the brane metric only. 



B. The gravitational field equations 



The final, but the most important constraints on the metric follow from the 5D Einstein 
equations. The independent components of the Einstein tensor of the metric given by Eq. ([2]) 
are given by 



V P^N^ PN^Q J ' \P PQj ' V PQ^ P^Q^ J 



P,y- 



N^y^yM^ 2P^r,rM^ '2P,y,yM^ Q,r,rM^ M 



2 



NQ^ PiV2 PQ2 iV2Q~ P^' ^ 

^''-ypM^MQj^^^^ypN^NQ)^''-^ — ^' 

f Nt 2Pt\ [Ny 2Py\ ^ Nty 2Pty 

- [m^MQj^^'-^yi^^m^J^^'^yMPQ^-^)^''^ 
fP, 2QA _ fP^ 2Q,N' \ fP^ _ 2Q,yN' \ 

\P PQJ \P^M^ PM^QJ \P^Q^ PQ3 J 



2 



M^y^yN' _ 2P,t,tN' '2P,y,yN^ _ Q,t,tN' _ N_ 
MQ^ PM^ PQ^ M^Q P2' y ) 

_ M^r,y , N^yM^r , Q,rM,y 2P,r,y 2NJ^ 2Q^ 
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22 



/ P'^N 



+ 



PP,r 

PP 



P 



+ I — H — 



MN^Q J 

Q,yP' 



fP'^Nt PPt 



'OA 



M,y- 



PP 



Q,yP^\ , PQ,rP,r PQ,tP,t 



2i 



32 



Ny + 



ATS 



PPu 
M2 



M2Q 
PP 

+ 



PQ,yP,y 

32/ 



f PP, 



\NM^ NM^q) 



+ 



32/ 



+ 



Mg2 



Q 



tV^y _j_ Q ,r,r P Q,t,t 



G. 



44 



MA^3 PMN"^ J 
^P^_2N^^^ 



2 



M2g 



(31) 



+ 



2P, 



+ 



2N. 



\P^m PN^ 



yP^A^P PM^N 

2 OD ^2 q2 
P2 



mj \PN 



MN PM 
P 



M,y + 



+ 



,y 
P 



P,y + 



Mm M'^N PN"^ PM^ 
The energy momentum tensor of the bulk is specified by Eq. f[T4l) 

/ A/f2/, n n Q n \ 






(32) 



T, 



AB 









































Q^NYrP 












'2 sin^ I 




Ar2+y2Q2 




Ar2+Y2Q2 y 



(33) 



In a source free bulk, two intermediate constraints of the metric can be obtained by setting 
Goi and Go4 to zero, 



2^,^ , Q, 



PM MQ 



MN 



+ 



2Pf 



PM J^^"^ \NQ ' PQ 



2Pr 

PiV 



+ 



+ 



Q r 

iVQ 
2P„ 



Q,t 



'2P,r,t 


Q,r,t 


P 


Q 


Nt,y 


'^P,t,y 


N 


P 



0, 



0. 



By differentiating with respect to r both sides of Eq. ([5]) gives 



M^r + y,rM^y = e*a^$. 



(34) 
(35) 

(36) 



By ehminating M,, with the use of Eq. (!36|) from the junction condition Eq. (fT5|) we obtain 

kl{2p + 3p)MNQ Y^e'^dr^ 



M 



■.y 



+ 



(37) 
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^From Eq. (jT]) we obtain 

Pt = ar and = a — YrP.y, (38) 
while from Eq. fl2T]) it follows that 

N^t = aAf, N^r = aAf-r - Fr-A^,y, (39) 

and 

Q,t = aQ, Q,r = aQ,r - Y,,Q^y, (40) 

respectively. By substituting the above relations and Eqs. ([5]), ([6]), and ([7]) into Eq. f p^ . we 
obtain the following equation for the metric tensor coefficients in the bulk evaluated on the 
brane, 

hi/ u \ VP V n 1 V n . 

0. (41) 



Q 



M J ar Q 

By substituting Eq. fl38l) into Eq. f|T8|) gives the following partial differential equation for 
the metric coefficient P, 

_klp PNQ aY,Q^ 

III. PARTICULAR CASES OF EMBEDDING OF BRANE WORLD WORM- 
HOLES 

In the previous Section we have formulated the basic equations describing the embedding 
of a four- dimensional braneworld wormhole into a five dimensional bulk. In the present 
Section we consider some particular cases of embedding. 



A. Braneworlds with homogeneous energy - momentum tensor 

As a first case in our study we consider the special case in which the energy and the 
pressure of the matter on the brane does not depend on the radial coordinate r, and it is a 
function of time only. Then the local conservation Eq. (125|) requires p = — p or 9,.$ = on 
the brane. 
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1. p^-p 



In this case, since dr^ = 0, Eq. ( !36|) becomes a first order linear partial differential 
equation, 

M^r + y,rM^y = 0, (43) 

with the general solution given by 

M{t,r,y) = M{t,y-Y{r)), (44) 

where J\4(t, z) = Jli(t,y — Y{r)) is an arbitrary function of the arguments. Eq. (1571) becomes 

M ,(t, y - Y{r)) = (2p + 3p) , ."^^ (45) 

Unless 2p + Bp = 0, this equation implies that the following expression is a constant, which 
for simplicity is taken to be 1, 

^2 1. (46) 



In terms of the throat function h{r) the above equation can be written as 

a/1 — b{r)/r ^ ^ 

Independently of the choice of Y{r), the bulk metric diverges at the throat. Since Af and Q 
are arbitrary integration functions, we can choose Q = 1, thus obtaining 



= T^, and = 0, 



Using Eq. fH2]) . we find 



,y-—^v^-Kr)/r. (48) 



^From Eq. (HTl) we can obtain Q^y, while Eq. (|T71) determines A^^j,. The time and extra- 
dimensional evolution of the bulk geometry is described by Eq. ( 135|) . Therefore the boundary 
conditions for the off brane evolution of the wormhole are completely determined. 

2. p = —p 

The local matter conservation given by Eq. tells us that if p = —p then p t = 0. 
Hence Eq. (I25l) is automatically satisfied. The freedom for the choice of the bulk metric for 
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this type of matter is larger. For example, a realizable setup is given by the following brane 
geometry, 

= l-6(r)/r' ^'^^^ 
which means that the wormhole was created solely by the embedding, and 

A/" = Q = 1, (50) 

respectively. Eqs. ([21]) and (HUj) simplify the junction condition given by Eq. ( TT7|) into 



My and P„ can be written in a simple form as 



M.= {'^ + a.^)<=yi-K (52) 



^.-(^)Vi-!:. (53) 



Once M and are given, the bulk Einstein equation Eq. ( 1411) determines Q,?;- From Eq. ( 135|1 
the form of M can be determined. This completely solves the problem of the embedding of 
the brane world wormhole into the bulk. 



B. Wormholes on a static brane 



Static brane world wormholes have already been studied in 20|. We would like now to 
investigate what would be the embedding for this static wormhole supported by the off- 
set of a brane. The induced metric is static, corresponding to a = 1. Since p and p are 
independent of time, the junction conditions given by Eqs. ( fT5l) - ( ITSl) tell us that the bulk 
metric is also static. The local matter conservation given by Eq. f l24|) and Eqs. (IMj) - ( 135|) 
is also automatically satisfied. Eq. fl25|) becomes 

p'=-{p + pW. (54) 
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A static wormhole is compatible with Q = 1. Using the junction conditions we can work 
out the y derivative of the metric coefficient as follows 



N 



,y 



6 



+ 



klNp (1 + Y^JN 



2N3/2 



Y P 



kl 



In particular a vacuum brane solution with p = p = is given by 



■.y 



P. 



Y,rM^r 

Y^rN^r + Y^r^rN 
m + TYj ' 



iV2 



Using Eq. (I38l) we obtain 



Y. 



1 - 



h[r) 



(55) 
(56) 

(57) 

(58) 
(59) 

(60) 
(61) 



With the use of these equations in the (^44 and components of the bulk Einstein field 
equations, and limiting ourselves to the brane, one could obtain two second order non-linear 
ordinary differential equations for M and A^, which involve only the second derivatives with 
respect to r of the metric functions. Therefore a consistent solution for a wormhole for a 
static brane always exists. 



C. Wormholes on a scalar field filled brane 



The energy momentum tensor of a scalar field is given by 



(62) 



where V (</>) is the self-interaction potential of the field. In the space-time described by the 
metric given by Eqs. ©-(IT]), the brane energy momentum tensor of the scalar field is given 
by 



diag < e 



,2$ 



-2* ^ 12 



-0^ + V{<j)) 



2 2 
a r 



-2$ i2 



2 

a r 
' r — b 

2 2-2/1 

, a r sm u 



-2<I> i2 



-0^ - r(0) 



(63) 
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This form of the energy-momentum tensor allows us to introduce an effective energy density 
and pressure according to the definitions 

Pit, r) = i0(t)'e-2*('-'*) + Vim, Pit, r) = ^^(tf e-'"^^^'^ - V^(0(t)). (64) 

With this forms of the effective energy density and pressure for the scalar field the conserva- 
tion equation Eq. ( l25ll is automatically satisfied. If the scalar field is potential dominated, 
it satisfies an effective equation of state p = —p. Moreover, the evolution of the scalar field 
is governed by the Klein - Gordon equation 

^ -_d,{^g>^^d,4>)+V'i(j))=0, (65) 



V-9 

where y/—g is the square root of the determinant of the metric tensor. In the geometry of 
Eqs. (I5])-(I7]), the evolution equation of the scalar field becomes 

4>+ (^3^-^]<^ + ^'('^)^' = 0' (66) 

which is consistent with the substitution of Eqs. flMj) into the conservation equation Eq. . 
Since asymptotically the wormhole geometry tends to the global (cosmological) geometry of 
the spacetime, by taking the limit r -> oo in Eq. ( l66l) . we found that the evolution of the 
scalar field (p results from the global cosmological evolution, 

+ 3-0 + \/'(0) = 0. (67) 
a 

Therefore we obtain 

M,(t,r,y(r)) ^ V'i<j)) 
M(t,r,F(r)) 

This equation can be integrated by using the substitution x = M^, and its general solution 
can be written as 



[M\t,r,Yir)) -1] . (68) 



where 



Ait) = (70) 



and M-ir) is an arbitrary integration function. Using Eqs. fj35|) and f HT]) a consistent set of 
M, A^, P and Q can be obtained. Therefore, an inflating brane world wormhole that does 
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not perturb the local evolution of the cosmological scalar field does exist. However, instead 
of studying its details by solving the close system of bulk Einstein equation, in the following 
Sections we derive some important results on brane world wormholes by using the effective 
brane description of the system. 



IV. THE BRANE GEOMETRY OF STATIC WORMHOLES 



In the covariant formulation of the RSII brane world model, the effective Einstein equation 



on the 4D brane is given by 



13| 



G^lu = Svr ^T^i, + ^'S'^;^^ - £^lu, (71) 

where A is the brane tension, T^j^ is the energy momentum tensor of the matter on the brane, 
and 

Sfj.1/ = -^TT^iy — —T^aT^ H dfiu, (72) 

is a quadratic term in the energy momentum tensor, which follows from the junction condi- 
tions of the embedding of the brane to the bulk. T = is the trace of the energy-momentum 
tensor. On the other hand, Ef^u = CABCDn'-'n^ g^g^ is the projection of the bulk Weyl tensor 
Cabcd to the brane. The metric (7^^ on the bulk induce a metric gj^J on the brane by the 
map that embed the brane. The metric of a static wormhole on the brane is given by j2o| 

ds^ = -e^^Mrft^ + + r\d9' + sin^ ed''cf>), (73) 

1 — o(r) /r 

where $(r) and b{r) are called the redshift function and the form function, respectively. To 
obtain a wormhole solution, several properties need to be imposed, namely [l|: The throat 
is located at r = ro and 6 (rg) = tq. A flaring out condition of the throat is required, i. e., 
{b — b'r) /b"^ > 0, where the prime denotes the derivative with respect to r. At the throat 
this inequality reduces to b' (tq) < 1. The condition 1 — b/r > is also required. To be 
traversable, there must be no horizons present, which are identifled as the surfaces with 
e^* — i- 0. Therefore $(r) must be flnite everywhere. For this metric the Einstein tensor 
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G 



diag 



1 



^2 ' ^ 



2$' h. h 



6' 



2 • 2 

, r sm I 



6' 



2r 



2r^ 



(74) 



The field equations on tlie brane can be written as G^y = SvrT^^, where T^^ = T^^ — 
{\/%-K)e^y + (6/A) Sfj,y. Once the energy momentum tensor on the brane T^^ is known, we 
can obtain the projected Weyl tensor by e^i, = SttT^^* — G^^, where T^^^ = T^^, + {Q/X)S^y. 
Since e^y is traceless, taking the trace of the brane Einstein equation Eq. (ffTl) would give 
the constraint equation SvrTrT^^* = TiGj^n, = —R, which relates the wormhole redshift and 
form functions with the matter component. From the trace of the Einstein tensor we obtain 



X 26' 
r + 36 - 4r) + — . 



This will equal to 



SttT = Stt 



(p-3p)- — (p' + 3p2--(p_3p)' 



(75) 



(76) 



In an orthonormal reference frame the components of the projected Weyl tensor e^/) have 
the components 

ef,o = diag [e(r), a,.(r), at(r), at{r)] . (77) 

The components of the effective energy-momentum tensor have the form p'^^ = p (1 + p/2X) — 
e/Svr, pf =p{l+ p/X) + pV2A - ar/Svr and pf =p{l + p/X) + p^/2X - at/Svr, respectively. 
The NEC violation, p'^^ + pf^ < provides the following generic restriction. 



871 (p + p) (l + jj< e + (Jr. 



(78) 



Hence braneworld gravity provides a natural scenario for the existence of traversable worm- 
holes |20|. 



V. 4D STUDY OF THE INFLATING BRANEWORLD WORMHOLE MODEL 

Roman [lo|] has suggested that inflation might provide a natural mechanism for the 
enlargement of Planck size wormholes to a macroscopic size. In the following we consider 
the inflation of the wormholes in the framework of the brane world models. 
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A. Brane geometry and the energy-momentum tensor 



The metric of an inflating wormliole on the brane is given by 



+ r\de'' + sin" Od^' 



(79) 



1 — b{r)/r 

In Eq. ( 1791) we have introduced a dynamical redshift function, but we do not consider the 
form function to evolve with time. This is because a nonzero db{rQ)/dt would require an 
infinite energy density of the matter at the throat. With this metric, the components of the 
Einstein tensor are given by 



Gtt 



+ 



b' 



G 



2d$' 



tr 



Gr 



r — b 



2$', b. b 
+ 1-- -- 



(80) 

(81) 
(82) 



Gn 



-2<S> 



[2ad^ - 2ad - d^] + ($" + $'2) " + 



2r V r 



(83) 



G, 



sin^ 6 



[2ad<i> - 2ad - d^] + ($" + ^'^) {^-^ + 



r \ 2r 2 1 2r^ \ r 



(84) 



We assume that the inflation is driven by a homogenous scalar field 0, filling the Universe, 
and with energy- momentum tensor given by Eq. fl62l) . The total matter contribution 
Tfj,i, + {6/\)S^u on the brane to the energy- momentum tensor in Eq. (17T|) is given by 



"Tnat 



mat 



mat 



2* , 

e e 



2 

a r 
r — b 



2 ^ ^ 8A ^ 2A ^ 2A 



0]/ _ 

2A ~ 2A ' ' 



-I mat 



2 2 I -2<i>^^ , -4*30^ , -2$0^ 

a r e V + e — — + e — r , , 

2 8A 2A 2A ' ' 



Z 



mat 



"^""^'^ T-^ + ^ 8A+^ 2A-2A 



(85) 
(86) 
(87) 
(88) 
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B. The nonlocal projected Weyl tensor 



The components of the projection of the bulk Weyl tensor Cabcd on the brane can be 
obtained from Eqs. flHU]) - flMl) and Eqs. flHSl) - flSSl) . respectively. These equations give the 
explicit form of e^^, in terms of the wormhole's redshift and form functions, and of the scalar 
field, respectively. The tt, tr and rr components of e^u are given by 



o 2<J> / I T/ I -4<I> '/'^ I -2<I> 

= 87re I e y + ^ + e ^ + e y^ 







)- 









' e2* 6' 



+ 



a? r2 



2(2$' 



(89) 
(90) 



r-6 2 ^ 8A ^ 2A 2A 



r — b 



2$ade-2* - 2aae-2* - d^e'^* + 



2$' 

r 



r 



(91) 



Taking the trace of the brane Einstein equation Eq. (17T1) gives another constraint equation 
for the Ricci tensor, 

2" 



6e 



2<I> 


d 


- $- + 


(:) 




a 


a 






\ 


/2$" 


2$'2 






a2 



+ 87r 

$'2 



_4$ V . , 

e — + e 
A 



A 



A 



j.2^2 



(b'r + 36 - 4r) 



26' 

^2q2 



0. 



(92) 



C. The energy conditions 

Similarly to the case of the static redshift inflating wormhole with $(r, t) = $(r), in order 
for the wormhole to be asymptotically fiat the "flaring out condition" must be satisfied. The 
metric of constant time and 6 = 7t/2 slice in the space defined by Eq. fl79|l is 

a^{t)dr^ 



1 — b{r)/r 



+ a\t)r'dn' 



With the use of Eq. (l93ll we can formulate the flaring out condition as lOj 

6 — 6'r 



262 



> 0. 



(93) 



(94) 



The energy condition at the throat can be worked out by considering the scalar quantity 
lim^-^ro T^^W^W^ , where = T^y + f 5^1/ — ^ is the effective energy momentum tensor 
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on the brane. For a radial outgoing null vector = (e 



lim T^IJM/^M/'^ = 2 

r—^ro 



a \ a -a 
- - - + $- 
a / a a 



-2<I> 



■, 0, 0) we obtain 
b' - 1 



(95) 



Therefore the energy conditions also evolve with the cosmological expansion of the Universe. 



VI. BRANE EVOLUTION OF THE INFLATING WORMHOLE 



The evolution of the scale factor a{t) and of the scalar field 0(t) are governed by the global 
evolution of the Universe. For an expanding wormhole the asymptotic behaviors of the form 
function and of the redshift function are given by lim^._^oo b{r)/r = and lim,r^oo ^{r) = 0, 
respectively. We also impose the condition limr^aob' /r^ = 0, so that asymptotically the 
behavior of the metric of the wormhole is the same as the behavior of the metric of the 
global braneworld model. Taking these limits in the dynamical equations describing the 
scale factor and the scalar field evolution gives the basic equations describing the inflating 
wormhole on the brane as 

2 



Stt 



02 04 ^2y y2 

— + v + — + \ 

2 8A 2A 2A 



and 



+ 3-0 + l^'(0) = 0, 



(96) 



(97) 



respectively. Since the two variables in these equations a{t) and 0(t) are functions of t only, 
we expect that these two equations actually hold everywhere in the wormhole. 



The dynamics of a scalar field in the brane world cosmology were studied in ISj and 



281], respectively, by assuming that the scalar field is confined in the 4-dimensional world. 



As for the potential of the scalar field, several types of potentials were considered. It has 
been shown that when the energy density square term is the dominating term in the energy- 
momentum tensor, the behavior of the scalar field is very different from the conventional 
cosmology. In the following we will restrict the discussion of the inflating wormhole to the 
scalar field potential 

Vi<P) = \m^\ (98) 
which corresponds to the chaotic inflation model on the brane 
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A. Time evolution of the redshift function 



The dynamic of the redshift function of a braneworld wormhole is given by Eq. (!68l) . 
Using the expressions of the metric tensor components on the induced metric we obtain 



= T/'(0) (e^* - 1) . (99) 

For e^* 7^ 1, we can separate the variables of the above equation into 

$ V'{(j)) 



e2*-l 

With 



(100) 



(101) 



dx , / e^^ — 1 



q2x _ 1 y g2x 

and by using V{(l)) = m^cj)^ /2, we obtain for the redshift function the expression 

where $(r, to) is the initial redshift function, which must be finite everywhere. The asymp- 
totic behavior of the redshift function is given by limr^oo^{f,t) = 1, and this behavior is 
the same for all times. The time evolution of the redshift function at the throat is described 
by the equation 



e2*(^«'*) = — —, (103) 

1 - A(ro)e-'*o 

where A{r) = [l — e'^**^*"'*"-'] can be obtained from the initial conditions. By taking the 
limit of small A we obtain 

g2$(.o,io) = ^ 1 + ^(^^) ^ ^^^^^^ = ^(^^). (104) 

1 - A{ro) 

This equation shows that A[r) is actually the Newtonian potential, and therefore we assume 
that A{r) < 0. 

VII. A SIMPLE SOLUTION - THE ZERO REDSHIFT CASE 

If the initial condition of the wormhole is so that $(r, Iq) = 0, Eq. fl99l) implies that $ = 
for all times. Then we also obtain 

0. (105) 
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Since b{ro) = tq, the solution for the form function is b{r) = tq. For this case the non - zero 
components of the Einstein tensor are given by 



Gtt — 3 



2 

a r 



r — 



2 2 

a r 



-2- 



a a 



G^^ = a'r'sm\e) 



a" 












+ 


ro 




2a2r3 




a 


-- + 
a'' 




_ 

a 





(106) 
(107) 
(108) 
(109) 



The energy - momentum tensor of the scalar field on the brane world wormhole is obtained 



as 



1 



2 

a r 



T^, = diag -0^ + r(0), 

'2 r — ro 



2 2 

a r 



2'^' - yi^) 



a^r' sin' (6) 



(110) 
(111) 



Eqs. f l89|) and Eq. fl9T]) give the exact form of the projected Weyl tensor on the brane as 

ro ro ro 



= diag 0, 



(l _ n.)' 2r' 2r 



:il2) 



One can easily verify that g^^e^u = 0. In the local coordinate system, the projected Weyl 
tensor can be represented as e^^i, = diag (eo, Cr, et, et), with 



eo = 0, er 



ro ro 



We concluded that the Weyl contribution that supporting the wormhole decay as a"^. 



(113) 



VIII. CHAOTIC INFLATION IN THE BRANEWORLD MODEL 



According to 18|,|28[, in the case of the chaotic inflation driven by the scalar field potential 
V{(j)) = m'clP' /2, the expressions of a{t) and (j){t) can be approximated by two different sets 
of expressions, corresponding to the exponential inflationary period, and to the final stages 
of inflation, respectively. 
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A. The Universe during exponential inflation 



Let ti be the initial time of inflation, and let t/ to be the end time of the inflation. During 
inflation, the potential term dominates the kinetic energy term, so that V{(f)) >> 0^, and 
V » (j), respectively. Then Eqs. (1^ and (|97|) have the solutions 



a{t) 



3n) 

a I exp 



(ti - t) 



1/2 



(114) 
(115) 



where ti and ai are arbitrary integration constants, which can be obtained from the initial 
conditions, and m is the 5D Planck mass, respectively. The end of inflation is determined by 
the breaking down of the slow roll condition, V{(j)) ^ (f)"^ , ot V ^ (j). These two conditions flx 
the value of t/, which is given by the expression tj ^ ti — (2m) By using this expression 
for tf, we obtain the scale factor after the exponential inflation phase as 



af := a (tf) = ai exp 



3 ^ ' 



aiexp ( - — 



(116) 



Since ti << tf, if we adopt the e- folding number for the exponential inflation phase to be 
60, we obtain the following estimate for the ratio of the scale factors at the beginning and 
end of inflation, respectively, 

a{tf) 



a {ti 



exp 



1 

-12 + 



m^tl 



exp(60) 



;ii7) 



Therefore mti ^ V180.25 ^ 13.43, and mtf ^ 12.9. 



B. The oscillating scalar field phase 

After the end of the exponential inflation phase, the scalar fleld enters in an oscillating 
stage, and it keeps oscillating even after the radiation domination era. The oscillating 



solution is approximated by [l8 



A \ ■* sm{mt) 



mt^/'^ ' 



a2{t/t2f'\ 



(118) 
(119) 



where ^2 and 02 are arbitrary integration constants. We can choose ^2 to be the time for 
which the oscillating stage begins, and determine it from the initial conditions. For example, 
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if at the beginning of the oscillating stage the scalar field is one fourth of that at the end of 
the exponential inflation, then 

1/4A\^, 1 /A\^sin(mt2) 

iU) ^3^) 

Together with tf ^ ti — (2m) this equation gives the following algebraic equation for mt2, 

""^^ - sm\mt2) = (121) 
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By using Newton's method, we find that the zero of this equation next to is given by 
mt2 ~ 13.75. We will discuss later that the solution is insensitive to the way we choose to 
interpret the two stages. 

IX. THE EVOLUTION OF THE INFLATING WORMHOLE 

In the present Section we consider the evolution of the inflating brane world wormhole 
during the different phases of inflation 

A. Wormhole evolution during the exponential inflation phase 

With the use of Eqs. (11141) and Eq. (11151) . the integral in Eq. (11031) becomes 

/ ^= / m^{t-ti)dt= / {x-mti)dx. (122) 

Jti J t=ti J x=mti 

Thus, for example, a 60 e-folding exponential inflation gives the following solution, describing 
the time evolution of the redshift function of the inflating brane world wormhole 

e^*^'''*) = ^-^-2 -, for mt < mtf = 12.9. (123) 

1 - A{r) exp - mHit) ^ ' 

The time evolution of the energy condition at the throat during the inflationary phase can 
then be obtained from Eq. fl95|) and from the global evolution Eq. (I114p as 

lim Tf^W^W" = ^\\- A(ro)e'^-'"'*^*l - ^A(ro)(mt - mtO^e"^"™'*^* + 

r-->ro ^ 3 L J 3 

^l^^,{mt-mt^f for mt < 12.9. (124) 
afro 

In this phase, we see that the among of exotic matter is diluting as an exponential function 
of the time. Furthermore, there is a range of wormhole initial conditions such that the 
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wormhole could survive the exponential inflationary phase, that is, the energy condition 
on the effective matter remain after the exponential inflation. At the end of inflation, the 
energy condition is obtained by substituting mt = 12.9 into Eq. fll24p . 

zL-rn^ h' 1 

hm Tf,WWU=t, ^ ^ [1 - 9.3 X 10-^°A(ro)] + ^ < 0, (125) 

where we have also used the value of mti obtained earlier. For example, wormholes with 
shaping function constrained by the following equation could survive the exponential infla- 
tion, 

b' < -^^"-jfyi + 1. (126) 

B. Dynamical evolution of the wormhole in the oscillating scalar field phase 

To obtain the evolution of the wormhole near the end of the exponentially inflation phase, 
corresponding to t > 5 we substitute Eq. f lllSp and Eq. 01191) into the global evolution 
equation Eq. (11031) . The integral in Eq. (I103p becomes 

^ ' (127) 



t2 <p Jx=mt2 2x — tanx 
Recall that t2 is the time at which the oscillating stage begins. We can see that there is a 
singularity in the integral, no matter when we switch to the oscillating stage. The evolving 
wormhole will eventually face a collapsing phase. To numerically analyze the solution, we 
adopt the 60 e-folding and the initial condition given by Eq. (I120p . and we assume that 
the value of the redshift function in Eq. (I123p at the end of the exponential inflation is the 
initial value for this phase. Hence the evolution of the redshift function at the throat can 
be written as 



= —z r, for ^ < (128) 



1 



1 — A(ro)e~^°exp /^^^g j^4:Xtanxdx/ {2x — tan 



X] 



where Xc is the flrst singular point of the function 4xtana;/ [2a; — tana;] that is larger than 
13.75, i.e., the flrst root of the equation 2x — tan(x) = greater than 13.75. With Newton's 
method we obtain that Xc ~ 14.101725, which is e-folding dependent. Besides, the value 
e~^° is obtained from the exponent exp (m^t^/2 — 13.43mt) in Eq. fll23p at mt = 12.9. 

Based on the singular behavior at Xc of the integral in Eq. (I127p . we can derive the fate 
of a wormhole at the characteristic time x^ for different initial conditions A(ro). 
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a. A (ro) = 0. This initial condition describes the zero redshift wormhole. It is the 
only initial condition that makes $ to remain finite everywhere. However, due to the fact 
that the integral J^~^^r^^ Axtanxdx/ (2x — tanx) is unbounded, it follows that such an initial 
condition is very restrictive. It is because any small derivation from zero redshift would be 
enlarged. On the other hand, the energy condition for the effective matter on the brane is 
violated at the throat for all times. It is the Weyl contribution from the bulk supporting 
such violation. As we have seen in Section [Vllt this Weyl contribution decay as a"^. 

b. A (ro) < 0. This initial condition makes the redshift function e^* vanish, since the 
integral in Eq. f ll28p diverges. This case corresponds to the appearance of an event horizon 
at the throat of the wormhole, and it can be interpreted as a conversion of wormhole into a 



black hole 



22| . Once the wormhole is converted into a black hole, the metric Eq. fl7^ does 



not describe anymore the geometry. This result also shows that any dynamical wormhole 
that survives the exponential inflation phase would collapse to a black hole, immediately 
after the oscillating stage begins. To see how the energy condition evolves in this stage, we 
analyze the $ terms in Eq. ( I95l) . Since e^^* — ?■ oo, there must be one moment in which 
the 2<l>e~^* term diverges, and it is negative as $ decreases. Therefore the energy condition 
of the effective matter remain violated during the conversion, which may form a unique 
signature on this black hole. 



X. DISCUSSIONS AND FINAL REMARKS 



In the present paper we have considered the evolution of inflating wormholes in the 
braneworld model. As a flrst step in our study we have considered the conditions under 
which a braneworld wormhole can be embedded into a flve dimensional bulk. We have 
found that such an embedding is, at least in principle, always possible. As a general result, 
we have found that braneworld wormholes created in the early Universe would experience a 
transition to a black hole. If there is a large number of braneworld wormholes created in the 
early Universe, that would lead to the creation of a large number of black holes, and this 
would have signiflcant astrophysical implications for the global evolution of the Universe, 
including the large scale structure formation period. 

Soon after the big bang, due to the symmetry breaking, the Universe underwent a phase 
transition after which it was fllled up with a highly homogenous scalar fleld. However, 
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space-time inhomogeneities can still be generated in such a homogenous background. For 
example, the perturbations of the scalar field could generate the initial inhomogeneities that 
eventually would be responsible for structure formation. Another possibility for generat- 
ing some very early space-time inhomogeneities during inflation would be the presence of 
braneworld wormholes in the early Universe. The braneworld wormholes are solutions of the 
Einstein equations on the brane. In the braneworld scenario, the existence of these worm- 
holes does not necessarily leads to the violation of the weak energy conditions. Braneworld 
wormholes would evolve, and increase in size as the Universe expands. With the exception 
of the particular case in which the wormholes have exactly zero redshift everywhere, when 
the Universe switches from the exponential inflation to the oscillating stage, a very small 
initial value A (ro) of the redshift function would make the wormhole to collapse into a black 
hole. The transition wormhole - black hole is realized through the sudden appearance of an 
event horizon. Since these black holes formed at a very early stage in the evolution of the 
Universe, and since there are no lower limits on their size, the black holes created from a 
wormhole should be considered as primordial. Once a wormhole becomes a black hole, it 
will soon lose its co-expansion with the Universe, and it can gain mass through relatively 



35 



36|. 



slow accretion during the radiation dominated era of the Universe 

Suppose that the Universe was fllled with wormholes that are separated by a distance rg. 
We would like to estimate the throat radius, so that these wormholes do not signiflcantly 
change the energy content of the universe. In order for this condition to be satisfled, we 
compare the total energy of a wormhole with the energy of the cosmological background. 
From Eq. (IHOl) we know that if we consider the average over whole space, the effect of creating 
a wormhole is similar to adding a extra energy term to the Universe that decays as 
If we impose the condition that this extra energy, necessary to create the wormhole, is 
larger than the energy corresponding to a wormhole fllled with the cosmological background 
energy, we obtain the condition 

-^a^(fr > I pa^d^r, (129) 



ro ^ ^ JO 

where p is the energy of the cosmological background. To obtain an estimate of the wormhole 
throat after inflation we use the form function of the wormhole for the dust solution, b' = 
'yr^/r'^, which gives 

47r^ > p-7crl (130) 
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Since < 7 < 1, in this order of magnitude estimation we take 7 1. If we adopt 



a reheating temperature of Trh = 10^ GeV^ 33|, [SJ], we obtain for the density of the 
cosmological background the approximate value p ~ 0.6 x 10^^ eV^. By choosing the end 
of inflation at a value of the scale factor of the order of a ^ 10"^" we can obtain an 

estimate of tq as 

r-o < ~ 10' cm. (131) 

The mass limit corresponding to a Schwarzschild's black hole would be 

1 km 

Mo < 4^ Me ^ 6.6 x 10=^^ g, (132) 

3 km 

corresponding to around 30 solar masses. 

In the present paper we have also explored some specific physical properties and charac- 
teristics of the inflating wormhole geometries, and we have found the expressions describing 
the properties of the wormholes during the different phases of the inflation. In a future work 
we will consider in more details the possible astrophysical implications of our results. 

XI. ACKNOWLEDGMENTS 

We would like to thank to the anonymous referees for comments and suggestions that 
helped us to significantly improve the manuscript. The authors would like to thank Dr. 
Francisco S. N. Lobo for useful suggestions and discussions. The work described in this 
paper was supported by a grant from the Research Grants Council of the Hong Kong Special 
Administrative Region, China (Project No. HKU 701808P). 



[1] M. S. Morris and K. S. Thorne, Am. J. Phys. 56, 395 (1988). 
[2] M. Visser, Lorentzian wormholes: from Einstein to Hawking, AIP Press, (1995). 
[3] M. S. Morris, K. S. Thorne, and U. Yurtsever, Phys. Rev. Lett. 61, 1446 (1988) 
[4] V. P. Frolov and I. D. Novikov, Phys. Rev. D42, 1057, (1990) 

[5] J. P. S. Lemos, F. S. N. Lobo and S. Quinet de Oliveira, Phys. Rev. D68, 064004 (2003). 
F. S. N. Lobo, Phys. Rev. D71, 084011 (2005); F. S. N. Lobo, Phys. Rev. D71, 124022 (2005); 
P. K. F. Kuhfittig, Phys. Rev. D73, 084014 (2006); F. S. N. Lobo, Phys. Rev. D75, 064027 
(2007); F. S. N. Lobo, Phys. Rev. D73, 064028 (2006); C. G. Boehmer, T. Harko, F. S. N. 

28 



Lobo, Class. Quant. Grav. 25, 075016 (2008); F. S. N. Lobo, Class. Quant. Grav. 25, 175006 

(2008); T. Harko, Z. Kovacs, F. S. N. Lobo, Phys. Rev. D78, 084005 (2008); R. Garattini and 

F. S. N. Lobo, Phys. Lett. B671, 146 (2009); T. Harko, Z. Kovacs, F. S. N. Lobo, Phys. Rev. 

D79, 064001 (2009); F. S. N. Lobo and M. A. Olivcira, Phys. Rev. D80, 104012 (2009); N. 

Montelongo Garcia and F. S. N. Lobo, Phys. Rev. D82, 104018 (2010); M. Jamil, P. K. F. 

Kuhfittig, F. Rahaman, and S. A. Rakib, Eur. Phys. J. C 67, 513 (2010). 
[6] M. S. Morris, K. S. Thorne and U. Yurtsever, Phys. Rev. Lett. 61, 1446 (1988). 
[7] S. W. Hawking and G. F. R. Ellis, The Large Scale Structure of Spacetime (Cambridge 

Univeristy Press, London, 1973), pp. 88-96. 
[8] F. J. Tipler, Phys. Rev. D17, 2521 (1978); T. Roman, ibid. 33, 3526 (1986); A. Borde, Class. 

Quantum Grav. 4, 343 (1987) 
[9] C. G. Boehmer, T. Harko and F. S. N. Lobo, Phys. Rev. D76 084014 (2007). 
[10] T. A. Roman, Phys. Rev. D47, 1370 (1993). 
[11] L. Randall and R. Sundrum, Phys. Rev. Lett 83, 3370 (1999). 
[12] L. Randall and R. Sundrum, Phys. Rev. Lett 83, 4690 (1999). 

[13] M. Sasaki, T. Shiromizu and K. Maeda, Phys. Rev. D62, 024008 (2000); T. Shiromizu, K. 
Maeda and M. Sasaki, Phys. Rev. D62, 024012 (2000); K. Maeda, S. Mizuno and T. Torii, 
Phys. Rev. D68, 024033 (2003). 

[14] A. Campos and C. F. Sopuerta, Phys. Rev. D63, 104012 (2001); C.-M. Chen, T. Harko and 
M. K. Mak, Phys. Rev. D64, 044013 (2001); C.-M. Chen, T. Harko and M. K. Mak, Phys. 
Rev. D64, 124017 (2001); J. D. Barrow and R. Maartens, Phys. Lett. B532, 153 (2002); 
C.-M. Chen, T. Harko, W. F. Kao and M. K. Mak, Nucl. Phys. B64, 159 (2002); T. Harko 
and M. K. Mak, Class. Quantum Grav. 20, 407 (2003); C.-M. Chen, T. Harko, W. F. Kao 
and M. K. Mak, JCAP 0311, 005 (2003); T. Harko and M. K. Mak, Class. Quantum Grav. 
21, 1489 (2004); M. K. Mak and T. Harko, Phys. Rev. D 70, 024010 (2004); T. Harko and 
M. K. Mak, Phys. Rev. D69, 064020 (2004); A. N. Aliev and A. E. Gumrukcuoglu, Class. 
Quant. Grav. 21, 5081 (2004); M. Maziashvih, Phys. Lett. B627, 197 (2005); M. K. Mak and 
T. Harko, Phys. Rev. D71, 104022 (2005); T. Harko and K. S. Cheng, Astrophys. J. 636, 8 
(2006); L. A. Gergely, Phys. Rev. D74 024002, (2006); C. G. Bohmer and T. Harko, Class. 
Quantum Grav. 24, 3191 (2007); M. Heydari-Fard and H. R. Sepangi, Phys. Lett. B649, 
1 (2007); T. Harko and K. S. Cheng, Phys. Rev. D76, 044013 (2007); A. Viznyuk and Y. 

29 



Shtanov, Phys. Rev. D76, 064009 (2007); Z. Kovacs and L. A. Gergely, Phys. Rev. D77, 
024003 (2008); T. Harko and V. S. Sabau, Phys. Rev. D77, 104009 (2008); Z. Keresztes and 
L. A. Gergely, Ann. Physik 19, 249 (2010); Z. Keresztes and L. A. Gergely, Class. Quant. 
Grav. 27, 105009 (2010); L. A. Gergely, T. Harko, M. Dwornik, G. Kupi, and Z. Keresztes, 
to appear in MNRAS, iarXiv:1105.0159l (2011). 

[15] K. Umezu, K. Ichiki, T. kajino, G. J. Mathews, R. Nakamura, and M. Yahiro, Phys. Rev. 
D73, 063527 (2005). 

[16] P. S. Apostolopoulos and N. Tetradis, Phys. Rev. D71, 043506 (2005). 

[17] T. Harko, W. F. Choi, K. C. Wong and K. S. Cheng, JCAP 0806, 002 (2008). 

[18] S. Mizuno, K.-I. Maeda, and K. Yamamoto, Phys. Rev. D67, 023516 (2003). 

[19] K. C. Wong, K. S. Cheng and T. Harko, Eur. Phys. J. C68, 241 (2010). 

[20] F. S. N. Lobo, Phys. Rev. D75, 064027 (2007). 

[21] J. P. Leon, Mod. Phys. Lett. A16, 2291 (2001). 

[22] P. K. F. Kuhfittig, Schol. Res. Exch. 2008, 296158 (2008). 

[23] A. R. Liddle and A. J. Smith, Phys. Rev. D68, 061301 (2003). 

[24] P. Bowcock, C. Charmousis and R. Gregory, Class. Quantum. Grav. 17 4745, (2000). 
[25] S. Mukohyama, T. Shiromizu and K. Maeda, Phys. Rev. D62, 024028 (2000). 
[26] P. Binetruy, C. Deffayet and D. Langlois, Nucl. Phys. B. 615, 219 (2001). 
[27] W. Israel, Nuovo Cimento B44, 1 (1966). 

[28] R. Maartens, D. Wands, B. A. Bassett and I. P. C. Heard, Phys. Rev. D62, 041301 (2000). 

[29] Y. Himemoto and T. Tanaka, Phys. Rev. D67, 0844014 (2003). 

[30] J. H. Brodi and D. A. Easson, JCAP 0312, 004 (2003). 

[31] M. Sami, N. Dadhich and T. Shiromizu, Phys. Lett. B568, 118 (2003). 

[32] J. Polchinski, Phys. Rev. Lett 75, 4724 (1995). 

[33] R. H. Cyburt, J. Ellis, B. D. Fields and K. A. Olvie, Phys. Rev. D67, 103521 (2003). 
[34] J. R. Ellis, J. E. Kim and D. V. Nanopoulos, Phys. Lett. B145, 181 (1984). 
[35] B. J. Carr and S. W. Hawking, Mon. Not R. Astron. Soc. 168, 399 (1974). 
[36] B. J. Carr, T. Harada and H. Maeda, Class. Quantum. Grav. 27 (2010). 



30 



